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Abstract 

By performing time-of-fiight neutron scattering measurements on a large amount of single crystals 
of Sr2RuC>4, we studied detailed structure of the imaginary part of the dynamic spin susceptibility 
over a wide range of phase space. In the normal state at T = 5 K, strong incommensurate (IC) 
peaks were clearly observed at around Q c = (0.3, 0.3) up to at least Hlj = 80 meV. In addition, our 
data also show strong magnetic fluctuations that exist on the ridges connecting the IC peaks around 
the (ir, ir) point rather than around the F point. Our results are consistent with the semi-mean-field 
random phase approximation calculation for a two dimensional Fermi liquid with a characteristic 
energy of 5.0 meV. Furthermore, the IC fluctuations were observed even at room temperature. 

PACS numbers: Valid PACS appear here 
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Despite lots of theoretical and experimental studies done since the discovery of the super- 
conductivity below T c = 1.5 K-£, the nature of the superconductivity in ruthenate, Sr2Ru04, 
is still controversial. The electronic structure of the ruthenate has been well studied theo- 
retically-^ and experimentally-. The ti g electrons of Ru 4+ ions form three bands near the 
Fermi surface; d xz and d yz orbitals form quasi-one-dimensional a and ft sheets, while d xy 
forms a two-dimensional 7 sheets. Initially, it was believed that the a and ft sheets are re- 
sponsible for the magnetic response, and the 7 sheet for the superconducting mechanism^. 
The former sheets generate the incommensurate (IC) antiferromagnetic fluctuations, while 
the latter appears as the ferromagnetic fluctuations. 

The superconducting order parameter was proposed to be a p-wave state with a chiral 
p x ± ip y symmetry^. The origin of the p-wave superconductivity was thought to be the 
ferromagnetic fluctuations due to the 7 sheets. However, the imaginary part of the dynamic 
spin susceptibility, x"> measured by inelastic neutron scattering measurement has not ex- 
hibited any detectable ferromagnetic fluctuations in both the normal and superconducting 
states^ - —. Instead, in both states gapless IC spin fluctuations were observed up to 40 meV 
at Q c = (0.3, 0.3) which come from the Fermi surface nesting of the a and ft sheets^. 
Furthermore, the direct experimental tests to observe the expected broken time-reversal 
symmetry has been conflicting. Muon spin relaxation^, polarized neutron scattering exper- 
iment^ and Kerr effect^ measurements have been reported to support the p x ± ip y paring 
scenario, while scanning Hall bar and scanning superconducting quantum interference device 
measurements^ gave no evidence for the edge supercurrents that should be generated by 
such a broken time-reversal symmetry. 

Very recently, theoretical work based on the weakly coupled bands reported that the 
quasi-one-dimensional a and ft sheets rather than two-dimensional 7 sheet are responsible 
for the triplet superconductivity and also for the lack of the edge supercurrents^. Since the 
a and ft sheets generate strong magnetic fluctuations, it would be crucial to obtain detailed 
information about dynamical spin susceptibility over a wide range of the momentum (Q) and 
energy (hw) phase space. Many theoretical models^^— predict the IC spin fluctuations, 
but the fine structure of x" varies in different models; while the position of the main peaks 
in x" is quite stable with regard to the computational details, the exact shape of the two- 
dimensional ridges connecting these peaks depends very much on the exact degree of the a/ ft 
nesting, including the residual three-dimensional dispersion. As a result, in some papers the 
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FIG. 1: (color online). Contour map of the neutron scattering intensity as function of H and hw at 
T = 5 K. The incident neutron energies were (a) E{ = 101.6 meV and (b) 151.2 meV. <S(Q, Huj) was 
integrated along the K-direction from ±0.15 to ±0.55 and L-direction from —9 to 9, respectively. 
Arrows represent the IC peak positions. 

ridges are largely suppressed around T—^-^-, while in others with different parameters and 
an additional spin-orbit interaction they extend over the entire length of the Brillouin zone^. 
Previous neutron scattering experiments have been limited to obtain the detailed structure 
of x" i mainly due to small amounts of sample and also due to the triple-axis spectroscopy 
(TAS) that was utilized to cover only a limited area of the phase space. 

Here, we report our time-of-flight (TOF) neutron scattering experiment on ~ 30 g single 
crystals of Sr2Ru04 to show detailed structure of the spin fluctuations over a very wide 
range of Q — ftw phase space in the normal state of Sr2RuC>4. Our results clearly show 
that the strong IC peaks exist up to at least 80 meV and can be well explained by the 
phenomenological response function for a Fermi liquid. Our analysis yield the characteristic 
energy for the spin fluctuation, hw^Y of 5.0(2) meV; Tiujsy is close to k B T FL below which 
Sr 2 Ru04 exhibit two-dimensional Fermi liquid behavior in the resistivity and specific heat 
measurements^— . Our data also shows that the weaker ridge scattering exists around the 
(tt, 7t) point rather than the T point. Furthermore, unlike in the previous neutron study-^ 1 ^, 
the IC fluctuations were observed even at room temperature. 

Five single crystals of Sr 2 Ru0 4 with a total mass of 29.15 g were prepared by a floating- 
zone method^ 1 ^. They were co-aligned for our neutron scattering measurements performed 
at the TOF spectrometer 4SEASONS 3 ->2i located at J-PARC, Tokai in Japan. The crystals 
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were mounted in a way that the crystallographic c-axis was along the incident neutron 
beam, which allows one to probe the scattering in the (if K 0) plane for this quasi-two- 
dimensional system. The crystals were put into an aluminum sample can that was then 
attached to a closed-cycle displex refrigerator, and the measurements were done at 5 K and 
300 K. The incident neutron energies were selected by a Fermi chopper and a combination 
of two disk choppers^. Due to the unique data collection system^, for a given setup 
of the choppers, the scattering events from neutrons with different incident energies, E- u 
are collected simultaneously by the detectors. We selected two setups that collected the 
following sets of incident energies: one with E[ = 12.6, 21.6, 45.5, 151.2 meV and another 
with E x = 28.1, 48.3, 101.6, 336.6 meV. Here, we show the data obtained with three most 
relevant incident energies, E[ = 12.6 (for magnetic fluctuations), and 101.6 and 151.2 meV 
(for the overall view of the magnetic fluctuations and lattice vibrational modes). 

Figures Ufa) and [TJb) show the scattering intensity with energy transfer up to fkj = 
100 meV along the H direction. There are strong bands of scattering centered at around 
20 meV, 40 meV, 60 meV, and 85 meV. They are phonons that have already been reported 
by the previous neutron scattering study using a thermal triple-axis spectrometer (TAS}~ 
In addition, there are strong spin fluctuations at the IC positions centered at Q c = (0.3, 0.3) 
that exist up to fvjj of at least 80 meV (see the white arrows), which are consistent with a 
previous polarized inelastic neutron scattering experiment that reported the IC peaks up to 
40 meVA 

The magnetic fluctuations can be more clearly seen at low energies where no strong 
phonon modes exist. Figure E{a) shows the data up to hw = 10 meV. The strong IC peaks 
are centered at the characteristic wave vector transfer, Q c = (0.3,0.3). The peak position 
does not change with hu up to 80 meV (see Fig. [TJ). The full-width-of-the-half- maximum of 
the magnetic scattering along the if -direction does not change at low energies (see Fig. EJ^a)), 
suggesting a very high stiffness of the magnetic fluctuations. This will be discussed in detail 
later. 

In order to investigate the Q-dependence of the low energy magnetic fluctuations, we 
integrated S(Q, fvjj) over hw from 2.5 to 7 meV and L from —5 to 5 at T = 5 K, and the 
result as a function of Q = (if, K) is plotted in Fig. [2]^b). There exist several IC peaks with 
the characteristic wave vector of Q c = (0.3,0.3). Positions of the IC peaks are summarized 
in the inset of Fig. [3J The IC magnetic fluctuations were theoretically predicted as the 
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FIG. 2: (color online), (a) Contour map of the neutron scattering intensity at T = 5 K as function 
of hco and H. The intensity was integrated over K from 0.25 to 0.38 and L from —5 to 5. (b) 
Contour map at constant- ftw cut the neutron scattering intensity at T = 5 K. The energy window 
was from hu> = 2.5 to 7 meV, and the scattering was integrated over —5 < L < 5. The incident 
neutron energy was E\ = 12.6 meV in both panels. 

consequence of the Fermi surface nesting of the quasi-one- dimensional a and f3 sheets^^i 
and have been observed experimentally^— . On the other hand, these peaks showed no 
L- dep endencoii^ 1 ^. 

To address the aforementioned controversy over the location of the ridge scattering, we 
integrated the neutron scattering intensity of the low energy excitations shown in Fig. Mjo) 
along several different Q-directions as shown in the inset of Fig. [3j Figure |3] shows the 
results. Figures [3]^a) and^b) show the strong IC peaks at Q c at 5 K (see red circles). Using 
the combination of gaussians and the linear background, the peak positions are obtained 
as follows: H = -0.311(1), 0.282(1), and 0.660(1) with K w 0.33, and K = -0.327(5), 
0.292(2), and 0.686(5) with H m 0.33. As shown as black triangles, the IC fluctuations 
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are broad but they survive even at 300 K. In addition, there are some shoulder peaks 
shown by the blue gaussians at H = -0.179(1), 0.175(3), and 0.720(4) with K w 0.33, 
and at K = -0.260(2), 0.283(1), and 0.781(9) with H « 0.33. These two features have 
already been observed in the previous neutron scattering studyii. For the possible ridge 
scattering, we plot the data across the ridges around (tt, tt) (see the inset for the directions, 
and Figs. E^c) and|3^d) for the data) and around the T point (Fig. EJe)). In contrast to the 
previous studyii, our TOF data show strong ridge scattering around (tt, tt) (at H = 0.326(1) 
with K fa 0.5 and K = 0.327(1) with H 0.5 in Figs. Ws) and[3^d), respectively) rather 
than around the Y point (Fig. EUe)). This can be seen in Figs. Ela) and[^b) as well; there 
are intensities between the peaks at H (or K) = 0.3 and 0.7 represented by green gaussian, 
but no additional intensity to the background between H (or K) = —0.3 and 0.3 except for 
the shoulder intensities. It is an unresolved question how the previous TAS experiment has 
yielded the different ridge scattering^. 

To study ^-dependence of the IC magnetic fluctuations at 5 K, we obtained the constant- 
fauu cuts at several different energy transfers as a function of H or K. Figure H] shows the 
results. Figures H(a) andHJ^b) show the low energy IC peaks can be easily extracted from the 
data shown in Fig. 2(a). Figures H]^c-h) show the IC peaks at high energies where phonons 
are weak and the magnetic scattering can be extracted from the data shown in Figs. 1(a) 
and (b). Our data clearly show that IC magnetic fluctuations in Sr 2 Ru04 survive all the 
way up to ~ 80 meV. The magnetic scattering intensity is related to the imaginary part of 
the dynamic susceptibility, x"(Q, hw), using the fluctuation dissipation theorem, 



The magnetic form factor of S^RuCi, F(Q), was obtained from reference^. For a quantita- 
tive analysis, we fit the data to the general form of the phenomenological response function 
used to describe a Fermi liquid system^; 



where x$> K o> ^sf, and Q c are parameters for the peak intensity, the sharpness of the 
peak, the characteristic energy of the spin fluctuations, and the IC peak position^ 1 ^, which 
are all independent of hw. We fitted the H- and independences at each hcu positions 
simultaneously to Eqs. (0Q) and © with a linear background. Q c was obtained to be Q c = 



S(Q,Hu,) = \F(Q)\ 2 



X"(Q,M 



(1) 



1 — exp (—hu I 'k^T) 
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(0.301(3), 0.303(4)) from fitting the 3 meV data (Figs. Ha) and Kb)), and fixed for all 
other energies. As shown by the solid lines in Fig. HI Eq. ([2]) describes our data well 
over the entire energy range. The optimum parameters obtained from the best fit are 
Xs = 10.6(3), k = 0.053(2) r.l.u. (= 0.086(4) A" 1 ) and frw SF = 5.0(2) meV. The fact 
that hwsF/kB = 58(3) K is close to Tfl may indicate that the two-dimensional IC spin 
fluctuations may be responsible for the Fermi liquid behavior. 

Our inelastic neutron scattering measurements on Sr2Ru04 over a very wide range of 
the Q — hu phase space revealed the three components of the magnetic fluctuations; the 
strong IC spin fluctuations at 5 K centered at Q c = (0.301(3), 0.303(4)) that extend up to at 
least 80 meV with the characteristic energy of hw^F = 5.0(2) meV, and the weaker so-called 
shoulder scattering and the ridge scattering at low energy transfer. Our data clearly shows 
that the ridge scattering is strong around the (it, tt) rather than around the T points. We 
have also shown that the IC fluctuations survive up to 300 K that is much higher than 
T c . Together with the previous observation that the magnetic fluctuations do not change 
when the system enters the superconducting phase below T^, our results suggest that the 
magnetic and superconducting properties of Sr 2 Ru04 are not connected. 

We thank Y.-B. Kim and I. Mazin for helpful discussion, and Y. Qiu for his help with 
analyzing the data. Works at the University of Virginia were supported by the US NSF 
under Agreement No. DMR-0903977. 
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FIG. 3: (color online). Several constant-Zio; cuts of the neutron scattering intensities at T = 5 K (red 
circles) and 300 K (black triangles) along different directions in the reciprocal space as described by 
arrows in inset: along (a) (H 0.33 0), (b) (H 0.5 0), (c) (0.33 K 0), (d) (0.5 K 0), and (e) (0 K 0). 
The energy window was from tvjj = 2.5 to 7 meV, and L was integrated from —5 to 5. The incident 
neutron energy was E\ = 12.6 meV. Striped areas (green) represent the ridge intensities, and solid 
areas (blue) represent the shoulder peaks. Solid lines are fitting results using the combination of 
gaussians and the linear background. The inset describes the schematic view of the (H K 0) plane 
in the reciprocal space. Big solid circles (black) represent the nuclear Bragg reflections. Small solid 
circles (blue) represent the IC magnetic peaks, thin lines (blue) show the positions of the shoulder 
peaks, and dashed lines (green) describe the ridge intensities between IC peak o 11 ' 14 ' 21 . Five thin 
arrows (red) show the directions (a, b, c, d, and e) in the reciprocal plane for each constant-foj 
cuts shown in Fig. [3l 
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FIG. 4: (color online), if- and ^-dependences of the neutron scattering intensity in the normal 
state at (a, b) Hui = 2, 3, 5, 7, (c, d) 32, (e, f) 55, and (g, h) 77.5 meV. H and K were integrated 
from 0.25 to 0.4, respectively, and L was integrated from —9 to 9. The incident neutron energies 
were (a, b) E\ = 12.6 meV, (c, d) 101.6 meV, and (e-h) 151.2 meV. Fitting results using Eqs. JT]) 
and (|2|) with a linear background at each hu are also described by solid lines. 
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